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We report measurements of the effects of a random vector potential generated by applying an in-plane mag- 
netic field to a graphene flake. Magnetic flux through the ripples cause orbital effects: phase-coherent weak 
localization is suppressed, while quasi-random Lorentz forces lead to anisotropic magnetoresistance. Distinct 
signatures of these two effects enable an independent estimation of the ripple amplitude and correlation length. 



Graphene is a one-atom-thick carbon sheet with unusual 
electronic properties due to its two-dimensional honeycomb 
lattice™. As an ultrathin membrane, graphene easily wrin- 
kles into the third dimension, with nanometer-scale ripples 
that have been observed in microscopic studies 3 6 . Strains as- 
sociated with rippling are expected to modify electronic trans- 
port by generating random scalar and vector potentials, which 
would affect transport by suppressing anti-localizatiorP and 
increasing the scattering rate^. Until this time, however, no 
transport measurements have directly probed graphene's rip- 
ples. 

Here, we report a magnetotransport measurement of 
graphene ripples using an in-plane magnetic field. In general, 
in-plane fields do not affect electron transport directly because 
electronic motion couples only to the component of magnetic 
field perpendicular to the graphene sheet. As illustrated in 
Fig. 1, however, a magnetic field aligned in the average plane 
of a graphene sheet (Sii) includes an inhomogeneous perpen- 
dicular component (5B±), which depends on the local slope 
of the graphene flake. Ripples allow the in-plane field to af- 
fect transport directly by converting it to an inhomogeneous 
out-of-plane magnetic field, a form of random vector potental 
(RVP). 

We observe two distinct effects of the Bu -controlled RVP 
that each depend on the ripples' rms height, Z, and correla- 
tion length, R. Random Aharonov-Bohm phases break time- 
reversal symmetry and suppress weak localizatiorPCS (WL) 
with an effective dephasing rate!^] proportional to B 2 Z 2 R. 

Random Lorentz forces^ lead to anisotropic momentum scat- 
tering by the in-plane field 14 and anisotropic magnetoresis- 
tance oc B 2 Z 2 /R. The effects are distinguished by measur- 
ing conductance in an additional uniform perpendicular field 
(B±) — an independent experimental knob that can be used to 
suppress time-reversal symmetry. 

Figure lb shows a schematic of the graphene flakes and 
measurement set-up. Two flakes (called A and B) were pre- 
pared on SiC>2, electrically contacted, then cooled to an elec- 
tron temperature of 40 mK. A two-magnet system provided 
independent control over the uniform magnetic field compo- 
nents Bii and B±. The carrier density n of the graphene was 
controlled capacitively using the Si back-gate (not shown in 
Fig. lb), and the electrical conductance G(n) showed a typical 
density-dependence (Fig. lc) with mesoscopic fluctuations. 
The conductance fluctuations of flakes A and B showed spin- 
split character when an in-plane field was applied, as reported 
in our earlier wortH Here, we instead report the magnetic 
field dependence of the mean conductance (G)„, averaging 
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FIG. 1: (a) Simulation of a rippled graphene sheet with a correlation 
length, R, ten times its rms height, Z. The uniform in-plane field 
B\\ , applied to the rippled topography of graphene, leads to a ran- 
dom surface-normal field SB±. (b) Schematics of graphene devices 
A (left) and B (right), showing orientation of applied fields Bu , B±, 
and electrical measurement setups (two-probe for A, multi-probe for 
B). Unused/broken electrodes are indicated by dashed edges. Scale 
bars are 5 /im. (c) Conductance G(n) for 5=0, flake B. (d) Weak- 
localization Bx-magnetoconductance of flake B at low density, mea- 
sured at different in-plane fields: zero (+), —4 T (□), +4 T (■). 
Two fits to {TJ are shown as solid lines. Applying the in-plane field 
dulls the central dip and decreases overall conductance. These are 
attributed to dephasing and scattering by the random 5Bj_. 



out (over intervals in n) the conductance fluctuations, follow- 
ing the procedure of Ref . |6) 

Figure Id shows typical (G) n (B±) curves for Bn = 
as well as B\\ = 4. We first discuss the conventional case, 
B\\ = 0. The conductance can be separated into Drude 
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and WL contributions, (G) n (B±) = ^[30 ± 9wl(B±)], 
with device aspect ratio ^r. The Drude conductivity, 3d = 
2e 2 vT m y/ir\n\/h, depends on the momentum scattering time 
r m and Fermi speed v « 10 6 m/s and does not change sig- 
nificantly with B± for the fields used in this experiment. The 
_Bj^-magnetoconductivity, ^[(G} n (B±) — (G)„(0)], is thus 
determined entirely by the WL component A(?wl(5_l)- 

Weak localization is a phase-coherent back-scattering ef- 
fect originating from the interference of a closed path with 
its reversed counterpart. Paths that contribute to WL are lim- 
ited in size by accumulated phases that are not symmetric un- 
der time (and path) reversal, coming from inelastic scatter- 
ing, magnetic fields and some types of disorder. Graphene's 
chiral charge carriers are intrinsically anti-localized, so 3wl 
is positive as long as chiral memory is maintained along a 
path. When chiral memory is lost due to intervalley scatter- 
ing, however, conventional weak localization decreases c/wl- 
Taking these effects together, a non-monotonic magnetocon- 
ductance curve is typically observed: large B± probes only 
short paths where chiral memory is retained, giving gwL > 0, 
whereas small B± also includes the contribution of longer 
paths for which chiral memory is lost, giving a conductance 
dip at B± = 0. 

For 2D diffusive graphene, the WL contribution has been 
calculated to beM^ 
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where F(z) = ln(z) +ip(~ + |) for digamma function i/j(x). 
Equation ([1} depends on four rates {r^ , t7 , Tj - , t^ 1 } 
characterizing different mechanisms that suppress WL. The 
diffusive accumulation of Aharonov-Bohm phases from uni- 
form B±_ gives Tg 1 = 2v 2 r m eBj_/h. The dephasing (t^ 1 ), 
inter- valley scattering (r" 1 ), and intra-valley scattering (r" 1 ) 
rates each originate from scattering processes that each break 
different time-reversal symmetries!^ 

The WL scattering rates were extracted from measured 
(G) n (B±) curves by fitting to Eq.[T] and are listed in Table]!] 
for the Bn = case. Values of r m (used to scale 1 ), as com- 
puted from 3d ~ x (G) n , were the primary systematic error 
since W- was difficult to determine (see Methods section). The 
longest time-scale, r^, may have been saturated by the device 
dimensions at high doping, with = vU T m T^j2 ~ 2/zm. 
The other characteristic lengths L,; ~ 600 nm, ~ 100 nm, 
and vr m < 100 nm, were not influenced by the sample geom- 
etry. 

Adding an in-plane field, Bu = 4T, changed (G) n (B±) in 
two distinct ways (Fig. Id). First, the dephasing rate was in- 
creased, visible as a suppression of the WL dip at small B±. 
Second, the overall (Drude) conductance was reduced, caus- 
ing an overall downward shift in the conductance at large B±. 
These effects can both be attributed to the ripple-induced RVP, 
and are discussed in order. 
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Low density 


Electron 


Hole 
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1 nil l „'2. 

11) /cm 


-13...-5 


-2.. .2 


5...13 


-5...-3 


W/L 




2.4 ± 0.8 


2.0 ± 0.7 


1.6 ±0.6 


0.7 ±0.3 


' m 


10 12 /s 


15 ±5 


20 ± 10 


15 ±5 


11 ± 5 




10 9 /s 


11 ± 1 


35 ±8 


11 ± 1 


11 ± 2 


$ 


10 9 /s 


70 ±50 


170 ± 70 


120 ± 80 


20 ± 10 


r- 1 


10 12 /s 


5.3 ±0.4 


2.7 ±0.5 


2.1 ±0.4 4.0 ±0.3 



TABLE I: Rates extracted from (G)„(Bj_, By = 0) at 40 mK (e.j 
Fig. Id); (G) n was averaged over the specified density ranges. 
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FIG. 2: (a,b) Bj_-magnetoconductance for small B±, at various val- 
ues of By. Fits toEq.[T] were computed assuming t.;, t», and r m are 
independent of By (as in Table while t^ 1 was a free parameter, 
(a) and (b) correspond to the low density (Fig. Id) and hole-doped 
regions, respectively, (c) For both flakes, extracted values of t^ 1 
increase in proportion with Bjj, as predicted from Eq. |2] Error bars 
plotted here do not include the larger systematic error from W- . 



The additional dephasing effect of an in-plane field due 
to Gaussian-correlated ripples was calculated in Ref. [T2] 
Whereas a uniform B± affects WL through the diffusive 
rate r^ 1 , the RVP affects WL as a micro-scattering rate 
(t^ 1 ) since the ripples are uncorrected beyond short dis- 
tances (R <C vr m ~ 100 nm): 

t- 1 ^t- 1 + V^(e 2 /h 2 )vZ 2 RBl (2) 
Eq. ([TJ was fit to multiple (G) n (B±) curves at finite By, al- 
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lowing only 1 to change from the B\\ — fits (Fig. 2a,b), 
in order to extract the Bn effect. Figure 2c confirms the 
Ar^ 1 (S||) cx B 2 dependence in Eq. |2j with a density- 
independent Z 2 R = 1.7 ± 0.5 nm 3 extracted for flakes A 
and B (Fig. 2); the uncertainty is dominated by uncertainty 
in Although the WL rate t^ 1 is commonly associated 
with inelastic scattering and loss of phase information, in this 
case it is enhanced by elastic scattering from the RVP which 
only scrambles phase information deterministically. The dis- 
tinction can be seen in conductance fluctuations, which are 
softened by decoherence but only scrambled by elastic scat- 
tering. As reported in our previous work on flakes A and B, 
fluctuations were scrambled by B\\ , and only decreased by | 
in variance due to broken spin symmetry.^ 

Besides suppressing localization, Bn caused an overall 
downward shift in the Bj^-magnetoconductance trace away 
from zero field (Fig. Id). This shift indicates a change in 
the Drude conductivity; it was isolated from the dephasing 
effect by examining the £?j|-magnetoconductance for values 
of\B±\ > 50 mT (Fig. 3), where the total WL correctiorP 
<?wl is essentially unaffected by the changes in . The ef- 
fect of RVP on Drude conductivity can understood as a de- 
crease in r m due to scattering by the Lorentz forces from the 
RVP (Fig. 3a). A similar random-field resistivity has been ob- 
served in 2D systems subject to RVPs originati ng from nearby 
magnetic particles or superconducting vortices^- 3 * 14 ! 17 ! 

The expected magnetoresistivity Ap(B\\) = 
^-[l/G(Bn) — 1/G(0)] for Gaussian ripples due to Lorentz 
forces from the RVP can be calculated by a Boltzmann 
approach, assuming kpR ^> 1 (high doping): 
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where 9 is the angle of the current flow relative to Bn . Equa- 
tion [3] is derived in the Appendix, and is most conveniently 
measured as a change in resistivity (rather than conductance). 
The threefold anisotropy predicted by Eq. Q was observed 
in a GaAs 2DEG when in-plane field lines were rippled by 
nearby ferromagnets 14 . In that case, as with rippled graphene, 
each individual magnetic ripple includes equal parts positive 
and negative magnetic flux, with a y-oriented zero-field chan- 
nel (for Bn along x; see Fig. 3a). In the strong-field limit, 
such channels might form magnetic waveguide d 18 * 19 ! but we 
observed a weaker form: the i -component of the velocity of a 
diffusive ensemble is randomized faster than the {/-component 
because {/-moving particles are deflected more (Fig. 3a)P^ 

The density dependence, Ap(n) cx |n|" 3 / 2 , predicted from 
Eq. [3] can be seen in the experimental data (Fig. 3b), on top 
of phase-coherent conductance (resistance) fluctuations due 
to the in-plane field 15 . The resistivity saturated for |n| < 
10 12 cm" 2 , perhaps due to the breakdown of classical scat- 
tering when kpR < 1. The fluctuations in Fig. 3b were 
minimized by measuring at 4K, and averaging the quantity 



%| 3/2 Ap(n, B\\) over n = (-3.5 . . . -1) x 10 12 cm" 2 , al- 
lowing fits to Eq.|3](Fig. 3c,d). It was confirmed that the mag- 
nitude of the effect did not change from 4K down to 40mK, 




FIG. 3: Anisotropic in-plane magnetoresistivity at B± > 50 mT 
in flake A at 4K. (a) Upper: a simulated symmetric bump in the 
graphene sheet, with uniform field Bu applied in the x direction. 
Lower: The resulting surface-normal field, SB±, is antisymmetric 
(positive on the right). Simulated trajectories show how an elec- 
tron's i-velocity is randomized more quickly than its y-velocity. (b) 
Density-dependence of p(8 T) - p(0 T) at B± = 50 mT, 6 = 20°. 
Predicted large-n behaviour in ^ for Z 2 /R — 0.15 nm shown as 
dashed curve, (c) Averages (Ii\n\ 3 ^ 2 Ap(n, B\\ ))„ show the depen- 
dence on the magnitude and direction of Bu . The current path (de- 
picted in inset) was measured at in-plane field orientations 9 m 20° 
(o) and 9 ~ 70° (+)• (d) Measurements on a different pair of elec- 
trodes (•¥ « 1.6) confirm that the anisotropy depends on the anglu- 
lar difference, 0, between field and current. 



though the amplitude of the fluctuations increased at low tem- 
perature as expected. 

Flake A was measured with two current paths along 9 w 
20° and 70°. The device was then re-cooled in a 90° -rotated 
orientation, to change 9 — > 9 + 90°. Fits of magnetoresistance 
curves to Eq.[3]gave a range Z 2 /R ~ 0.05-0.2 nm for flake A 
(Fig. 3cd). The measured anisotropy A / o(70 )/Ap(20°) was 
approximately 0.13 ± 0.01 for one current path (Fig. 3c) and 
0.26 ± 0.03 for the other (Fig. 3d), whereas Eq. ([3} predicts 
0.44. In the single measurement of flake B, Z 2 / R ss 0.02- 
0.04 nm. 

Using the value Z 2 R 1.7 nm 3 from the analysis in 
Fig. 2 and the range of values of Z 2 /R reported above, 
Z = 0.6 ± 0.1 nm and R = 4 ± 2 nm are extracted for 
Gaussian-correlated rippling of Flake A (the spread in Z 2 /R 
is incorporated into uncertainties for Z and R). The values 
for Z and R from this work can be compared to values ob- 
tained from scanning probe measurements on these flakes, 
and to values reported in the literature for graphene flakes 
on Si02- After the measurements described above, room 
temperature atomic force microscope (AFM) measurements 
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were performed on flake A using an Asylum MFP3D-SA, 
after annealing the flake at 400°C in a low pressure N2/H2 
gas mixture to remove resist residues 4 . These measurements 
gave Z = 0.13 ± 0.02 nm and R = 10 ± 5 nm; limi- 
tations of vibration and drift prevented more accurate mea- 
surements. AFM measurements on similar flakes in Ref. 4 
gave Z = 0.19 nm and R — 32 nm, whereas scanning 
tunnelling microscope (STM) measurements in Ref. |5] gave 
Z = 0.35 nm, R w 5 nrrPS. A consistent discrepancy is 
noted between scanning probe measurements and topographic 
parameters extracted from transport: the observed magneto- 
resistance in flakes A and B is nearly a factor of 100 larger 
than would be expected from our own AFM measurements on 
flake A. Although the right-hand side of Eq. ([3]) is enhanced 
by a factor between one and two for realistic non-Gaussian 
correlations in ripples (see Appendix), this is far too small to 
explain the measured magneto-resistance. 

The in-plane field couples to spins as well as to the mo- 
tion of charges, leading to the possibility of spin-related ef- 
fects on transport. 15 Three different types of spin-related ef- 
fects are considered. First, Zeeman splitting of the Fermi 
level implies altered populations of spin-up and spin-down 
electrons 15 that screen charged impurities less efficiently, 
leading to a magnetoresistance 2 -^ p(Bn)/p(0) sa (2 x 
10 5 cm~ 2 T~ 2 )i? 2 /|n| for densities larger than the impurity 
broadening. Quantitatively, however, this effect would be too 
small to be observed in our measurements. Second, spin- 
flip scattering off magnetic impurities can lead to decoher- 
ence, adding to . This effect would be field dependent, 

as B > Bi will freeze impurities into their ground state 22 , 
disabling spin-flip dephasing at B t — k^T/g*^ ~ 100 mT. 
Such an effect would show up as a peak in t7 at zero field 
in data such as Fig. 2c. Based on the apparent absence of 
this peak, impurity-induced spin flips must be uncommon for 
electrons traversing our devices. Finally, magnetic impurities 
may also generate an RVP by their localized magnetic fields, 
but the strength of this RVP (and hence its contribution to r^) 
would not change when the impurities align to B11 . 

Finally, we turn to an important analogy that can be drawn 
between the effects of an in-plane field and those of strain 
due to ripples. Since both the in-plane field and ripple strain 
generate random vect or po tentials that are directly correlated 
with ripple topography-^!, similar effects can be expected. In 
particular, strain is commonly believed to suppress weak anti- 
localizatiorPSSD but there is widespread disagreement about 
how to estimate the magnitude of the effect. We argue that the 
suppression occurs through a short-range dephasing process 
much like that in Eq.[2] and that strain-induced dephasing can 
fully explain the observed suppression of anti-localization. 

Strain RVP affects each valley oppositely with a magnitude 
that depends on R and Z™ corresponding to a fictitious valley- 
dependent in-plane field of size HZ/(eaoR 2 ), where the lat- 
tice constant is a = 0.14 nm. The resulting valley-dependent 
dephasing affects t~V 6 and by analogy with Eq. 2 we expect 
t^ 1 1=3 vZ 4 / (a 2 i? 3 ). Using this expression, the ripple dimen- 
sions extracted from our in-plane field measurements or from 
STM measurements 5 may fully explain the large r^ 1 ~ 1- 
10 ps _1 observed in most graphene WL magnetoresistance 



experiments (see Tableland Refs. 161231241) . 

This result contrasts with previous estimates of the intra- 
valley effect of ripple strain, which have assumed that the 
strain-induced effective magnetic field is truly random, with 
no requirem ent fo r flux compensation over multiple correla- 
tion lengthsP^HH That assumption would imply long range 
correlations in the RVP. The dephasing measurements in the 
first half of this paper show that ripples and resulting RVP cor- 
relations are instead short-range, as expected for adhesion to 
a polished wafer, and the analysis for in Eq. 2 should apply 
also to . 

In conclusion, transport measurements of graphene flakes 
in an in-plane magnetic field showed effects due to the mag- 
netic flux threaded through the ripples. The use of an auxil- 
iary out-of-plane field allowed two distinct effects to be sepa- 
rated: weak localization suppression (by dephasing) and over- 
all anisotropic magnetoresistance (by Lorentz-force scatter- 
ing). Besides allowing a determination of the ripples' typical 
height and length scale, these measurements provide insight 
as to how other short-range random vector potentials (such as 
that due to ripple strain) might affect transport in graphene. 



Experimental methods 

Silicon wafers with a ^300 nm wet thermal oxide were 
thinned to ^260 nm oxide thickness by a CF4/O2 plasma, 
then cleaned using a standard SC-l/SC-2 process. Within 
an hour of cleaning, flakes of graphene were deposited us- 
ing the mechanical exfoliation technique 2 then located in an 
optical microscope. The Cr/Au electrodes were deposited 
in an e-beam lithography process using PMMA resist. Im- 
mediately before the cooldowns for electrical measurements, 
devices were baked on a 125°C hot plate. Graphene flakes 
were confirmed to be single layers by quantum Hall effect 
measurements 15 ; the backgate capacitance ne/(Vc — Vq) = 
8.0±0.1 x 10 10 cm _2 e/V was calculated from resistivity min- 
ima in this data. In each cooldown, the gate offset Vq was 
determined by requiring the conductance minimum to occur 
at n = 0. For Figs. 1 and 2, V = 1 V in flake aP and 
Vo = 11.5 V in flake B. The quality of flake A had decayed 
before the 4 K measurements of Fig. 3, so that Vq varied from 
13 to 23 V depending on cooldown and the current path. 

Two magnets provided fields up to 120 mT oriented out- 
of-plane, and 12 T oriented roughly in-plane. A ~1° mis- 
alignment of the in-plane axis was corrected by biasing the 
out-of-plane magnet: the reported B± values have been man- 
ually offset at each Bt\ so that the WL dip always occurs at 
B± = 0. Conductance was measured by lock-in techniques 
with a 10 nA current bias (Fig. lb). In flake B, conductance 
was measured in a four-probe geometry (Fig. lb), and we fur- 
ther reduced the effects of conductance fluctuations by aver- 
aging over two opposing sets of voltage probes with the same 
aspect ratio (Fig. lb). In flake A, two-terminal conductance 
was calculated after subtracting a 3.2 kQ contact resistance. 
Aspect ratios of flake A were computed by a 2D Laplace equa- 
tion solver, and aspect ratios of flake B were estimated from 
the device geometry. In either case, aspect ratios were com- 
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plicated by the invasive nature of the contacts 26 . The differ- 
ent aspect ratios for hole and electron doping (Table [TJ> were 
determined by requiring the conductivity to be n-symmetric: 



(§G)(n) = (±G)(-n). 



Appendix: B\\ anisotropic magnetoresistance 

We model diffusive particle motion by a probability distri- 
bution f(k), with evolution ^ [/(£)] = D[S + S h f](k), 
for D[S, f](k) = / $frS(k, k>){f{k') - f{k)). The zero- 
field scattering matrix is So, assumed to have the elastic 
and isotropic form So{k,k') — So(k,q)5(k' — k), where 
q = k' — k. Transport occurs by modes f x (9) = cos(0) 
and fy(6) — -4= sin(6>) which are eigenmodes of D[Sq, /,*] — 



~ T u /« f° r a ^ The zero-field transport time is then t xx — 
T yy = T m , with T X y = 0. 

The Si I -induced scattering matrix is calculated by Fermi 

golden rule as S\\(k,k?) = ^f\(k'\U{r)\k)\ 2 5{E kl - E k ). 
The potential is U(r) = —evA(r) ■ a for Dirac fermions 
with Pauli matrices a operating on chirality. Given Sn ori- 
ented along x, the vector potential A(r) = —B\\h(r)y 
for rippled graphene with out-of -plane displacements denoted 
by h{r). We consider just single-valley chiral plane waves 



\k) 



-ie h /2 



/2V 



lfc / 2 ] T under the assumption that 



U does not mix valleys, yielding 

S\\(k,k') = fi2 l c (q) S m 2 ^^6(k-k'), (A.l) 

for height correlator c(r) — (h(fo)h(fo + r)). This scat- 
tering depends on the sum of 9k' and 9k, and therefore may 
anisotropically break the degeneracy of the transport eigen- 
modes. As a perturbation of the evolution, this modifies the 
scattering rates as St/j 1 = f* 7r d9 k fi(9 k )D[S\\, fj]{0k), to 
first order. 

Here we examine only the case of isotropic ripples c(q) = 
c(q), for which p xy ,p yx remain zero. After some calculation, 
the anisotropy is found to be exactly threefold: Ap xx = |py, 



and Apyy = I /On, where pi\ (k) may be written as a real-space 
integral of the height correlator: 



P\\(k) 



nB 2 



dr rW(kr)c(r), 



(A.2) 



for W(z) — dip Jo(2z sin |) sin 2 |. We have confirmed 
the threefold anisotropy of relaxation by simulating a classical 
charge moving in the x-y plane with out-of-plane magnetic 
field B = — SiixrZ, for random isotropic h(r). Simulation 
and experiment in Ref.fT4lalso show Ap xx = 3Ap yy . 

W(z) is approximately constant for z < 1, so pn (k) 
plateaus at a resistivity ~ (ZRB^) 2 /h at low densities where 
k <C l/R- At higher density (k > 1/-R), the oscillatory 
W(z) may be integrated out via a Hankel transform. This 
yields the classical result pn oc fc -3 of Eq. pjwith an effective 
[Z 2 /RU - - / °° dr i d. c ( r )/V5F, whicTTis exactly Z 2 / R 
if the ripples are Gaussian. 

For realistic non-Gaussian correlation functions, and for 
k « l/R, Eq.|3] may be used for Z and R statistically de- 
fined as in Ref. 4 though with some error. For more accuracy, 
Eq. A.2 should be computed directly. Using the correlation 



function of Ref.|4] Eq. A.2 gives a magnetoresistance approxi- 
mately 80% larger than that expected from Eq|3] due to fractal 
scaling of ripples which causes \Z 2 / R} c s > Z 2 /R. For the 
fiw5 nm correlation function 20 of Ref. 5J k l/R for the 
range of n examined in Fig. 3, and Eq. A.2 gives the magne- 



toresistance a fc-dependence slightly weaker than the classical 
k~ 3 . Averaging over n in the method of Fig. 3 yields a result 
in agreement with Eq. [3] though strong disagreement arises 
for R < 4 nm as pn saturates to its maximum. 
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